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SUMMARY 


This  report  gives  a  general  description  of  a  digital  computer  program 
used  in  connection  with  the  study  of  injury  of  skin  exposed  to  thermal  energy. 
All  of  the  information  necessary  for  a  detailed  understanding  of  the  program 
is  included;  however,  the  material  is  presented  in  a  manner  such  that  the  novice 
in  the  field  of  computer  science  may  make  use  of  the  program  if  he  so  desires. 
For  this  reason  emphasis  is  placed  on  the  operating  instructions  for  the  pro¬ 
gram.  A  short  discussion  of  the  pertinent  theory  and  equations  as  they  apply 
to  the  human  skin  is  included  at  the  beginning  of  this  report. 
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INTRODUCTION 


This  report  contains  a  description  of  a  digital  computer  program 
that  can  be  used  to  evaluate  theoretical  equations  associated  with  the  time- 
temperature  histories  of  skin  exposed  to  various  levels  of  thermal  radiation 
and  to  predict  the  injury  due  to  such  exposures  (1).  The  program  is  written 
for  a  Control  Data  Corporation  3200  Computer  System  using  the  Fortran  3200 
language.  In  the  study  of  thermal  tissue  damage  it  is  of  interest  to  obtain 
the  time-temperature  history  at  some  depth  below  the  surface  of  the  skin  such 
as  at  the  dermis-epidermis  interface  and  also  to  obtain  the  time-temperature 
history  at  the  surface  of  the  skin.  For  this  reason  the  computer  program  in¬ 
corporates  the  feature  of  obtaining  the  time-temperature  history  at  depth  and 
at  the  surface. 

THEORY  OF  PROBLEM  AND  EQUATIONS 

The  time-temperature  history  of  skin  exposed  to  a  square-wave  pulse  of 
thermal  energy  is  characterized  by  a  temperature  rise  from  some  initial  tempera¬ 
ture  of  the  surrounding  environment  T  at  time  t=o  when  irradiation  at  a  flux 
of  magnitude,  Q,  begins.  The  temperature  continue:  to  rise  to  some  peak  tempera¬ 
ture  at  time  t  *  i,  the  time  at  which  the  radiation  ceases.  The  temperature  then 
drops,  rapidly  at  first,  then  more  slowly,  approaching  Tq  as  t  approaches  in¬ 
finity.  Figure  1  shows  a  typical  example  of  a  t ime-temperature  history  as  com¬ 
puted  for  a  specific  exposure. 

The  following  equation  suggested  by  Buettner  (2)  describes  the  desired 
t ime-temperature  history: 
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Figure  1.  Typical  Time  -  Temperature  History  of  Skin  Exposed  to  a  Square- 
Wave  Pulse  of  Thermal  Energy. 
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Eq.  1 


where  Tx  »  tissue  temperature  at  depth  x  below  the  surface  (*C) 

2 

Q  ■  effect  it  radiation  on  the  surface  of  skin  (Cal/cm  sec) 
k  *  thermal  conductivity  (Cal/cm  sec  #C) 
p  *  density  of  skin  (g/cm^) 

c  *  specific  heat  (Cal/g  °C) 

2  2 
a  =  k/pc  =  "temperature  diffusivity"  (cm  /sec) 

t  «=  time  (sec) 

t  «=  time  at  which  thermal  radiation  ceases  (exposure  time)  (sec) 

x  =  depth  below  the  surface  of  the  skin  (cm) 

.0, 


Tq  *  initial  surrounding  temperature  (  C) 
8 (U)  *  integral  of  the  probability  curve  ■ 
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Function 


Equation  1  can  be  derived  directly  from  the  general  differential  equation  for  heat 
conduction  in  one  dimension 
fiT 


2  fi2T 
"  4  6x^" 


fit  ’  fix'  Eq>  2 

assuming  a  constant  heat  flow  and  an  initial  isothermal  condition  (vertical  temp¬ 
erature  gradient  equals  zero)  and  the  heat  absorbed  at  the  surface  of  the  skin 
transferred  inward  by  conduction. 

Prediction  of  dermal  injury  resulting  from  exposure  to  thermal  radiation 
of  any  given  magnitude  and  duration  depends  entirely  upon  the  resultant  tiae- 
temperature  history.  Total  tissue  damage  done  during  any  gven  episode  aust 


3 


include  the  damage  done  during  cooling  after  the  radiation  ceases  as  well  as 
the  damage  done  during  heating.  Equations  3  and  4  express  damage  as  temporal 
integral  of  rates  of  tissue  injury  depending  upon  the  tissue  temperature,  and 


necrosis 


dfl/dt  -  damage  rate  at  given  temperature,  Tx 

dt  =  time  interval  for  which  given  temperature  prevailed  (sec.) 

t  =  time  at  which  thermal  radiation  ceases  (exposure  time) (sec.) 

tj  =  time  at  which  the  injurious  temperature  level  *.44°C)  is 

reached  (sec.) 

ti  =  time  at  which  temperature  falls  below  the  injurious  level 
(44°C)  (sec.) 

P  =  constant  of  integration 

AE  =  energy  of  inactivation 

R  ■  gas  constant 

«  tissue  temperature  at  depth  x  in  %  at  time  t 
The  first  term  on  the  right  hand  side  of  Eq.  4  is  the  damage  done  during  heat¬ 
ing,  hereafter  designated  the  second  term  is  the  damage  dore  during  cooling, 
hereafter  designated  the  sum  of  these  two  terms  is  the  total  damage  done, 
hereafter  designated  ilj.  and  is  equal  to  unity  at  the  point  of  complete  trans- 
epidermal  necrosis. 
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PKOCRAM  DESCRIPTION 


The  complete  program  can  he  roughly  broken  into  seven  parts.  Figure  2 
is  a  listing  of  the  Fortran  statements  of  the  cntir  program. 

The  first  part  of  the  program  reads  into  the  computer  the  necessary 
data  required  for  calculation  of  t imc-temperature  histories.  This  includes 
data  such  as  the  number  of  time-temperature  histones  to  he  computed,  various 
labels  used  on  graph  outputs,  constants  of  inteerition,  etc. 

The  second  part  of  the  program  computes  the  individual  t ime-temperature 
histories  as  arrays  of  time-temperature  points,  kach  time  point  is  stored  in  the 
array  TIMEyN)  and  the  associated  temperature  point  is  stored  in  the  array  T(N). 

In  addition,  an  array  A(N)  containing  the  square  root  of  values  of  thermal  con¬ 
ductivity  is  stored.  Each  of  the  three  arrays  has  a  maximum  dimension  of  200 
floating  point  values.  For  each  time-temperature  historv  to  be  confuted,  a 
data  card  is  read  in  containing  the  following  variables: 

Q  =  effective  radiation  on  the  s-irface  of  the  skin 
X  =  depth  below  the  surface  of  sbn  (x  =  o  at  the  surface) 

DT  =  time  interval  bt"ween  points  in  the  t  ire-temperature  history 
TAU  =  time  at  which  thermal  radiation  ceases  =  t  'exposure  time) 

T IMF C 1 )  -  initial  starting  point  m  time  at  which  first  temperature  pcint 
is  computed 

A1  1  square  root  of  value  of  thermal  conductivity  during  heating  phase 
A2  -  square  root  of  value  of  thermal  conductivity  used  to  compute 
first  temperature  point  during  the  cooling  phase. 
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Figure  2.  Fortran  Statement  Listing  of  Program. 
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Continuation  of  Figure  2 
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Continuation  of  Figure  2 
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It  is  necessary  to  use  some  value  other  than  zero  for  TIME(l)  in  order  that  the 
value  of  the  exponent  containing  TIME(N)  as  a  divisor  does  not  approach  infinity 
since  in  this  event  an  exponent  overflow  error  would  occur.  Hence  a  value  of 
TIME ( 1 )  very  nearly  equal  to  zero  should  be  chosen  so  that  the  sum  of  TIME(l) 
plus  some  integer  number  times  DT  is  very  nearly  equal  to  TAU.  Again  it  is  neces 
sary  that  the  sum  of  TIME ( 1 )  plus  some  integer  number  times  DT  does  not  exactly 
equal  TAU  in  order  to  prevent  errcr. 

In  solving  Ej.  1  on  the  computer,  the  following  approximation  is  used: 


pc=  1 

Eq.  5 

Hence 

a2=  k 

or 

a  =  A 

Eq.  6 

The  computation  of  the  time-temperature  history  is  divided  into  two  intervals. 

0  <  TIME (N)  <_  TAU  and  TIME(N)  >  TAU,  which  correspond  to  the  heating  phase  and 
cooling  phase  respectively.  Before  the  calculation  of  each  temperature  point,  a 
check  is  ■"ade  to  see  if  TIME(N)  <_  TAU.  If  this  is  the  case,  Eq.  1  is  solved  for 
the  first  two  terms  on  the  right  hand  side,  the  third  term  being  imaginary.  When 
TIME(N)  exceeds  TAU  in  value  (during  the  cooling  phase),  the  entire  equation  is 
solved.  During  the  heating  phase  a  constant  value  of  the  square  root  of  thermal 
conductivity,  Al,  is  used;  however,  during  the  exiling  phase  the  value  of  the 
square  root  of  thermal  conductivity,  A2,  is  decremented  by  0.0001  after  the  cal¬ 
culation  of  the  first  temperature  point  and  is  continually  decremented  for  each 
temperature  point  thereafter 

The  program  uses  an  approximation  formula  to  compute  the  value  of  the 
complementary  Error  Function,  1  -  6(U).  The  approximation  formula  for  the  Error 
Function  8(U)  is  given  by  Hidings  (6). 


? 


Since  the  complementary  Error  Function  is  equal  to  one  minus  the  Error  Function 
we  have,  Con.  0(U)  ■  1-0 (U) 

1 

Co».  e<U)  -  COMERF  =  .  [l.aiu.a2U2.a3l|3.a4U4.a5y5.il()Ui]l6  tR-  8 


where 


aj  -  0.0705230784 
a2  *  0.0422820123 
a3  *  0.0092705272 
a4  *  0.0001520143 
a5  =  0.0002765672 
a6  =  0.0000430638 


During  and  following  calculation  of  the  t ime-ter.perature  history  the 
following  three  specific  values  of  T(N)  are  stored  separately  for  later  use: 

TPK  =  peak  temperature  obtained 

TMINUP  =  value  of  the  temperature  closest  to  44°C  (injurious  temperature 
level)  during  the  heating  phase. 

TMINDN  *  last  value  of  temperature  in  the  T(N)  array. 

Also  the  number  of  elements  in  the  array  T(N)  between  TMIM'P  and  TMINDN  inclusivt 
is  stored  in  NO  for  later  use.  Figure  3  is  a  flow  chart  of  time-temperature 
history  computation. 

The  last  five  parts  of  the  program  are  essentially  connected  with  out¬ 
put  and  are  selected  by  means  of  sense  switches  located  on  the  computer  console. 
Thus  any,  all,  or  none  of  the  five  parts  can  be  selected  in  turn.  Ihe  operation 
of  any  one  of  the  parts  may  be  skipped  in  the  program  ly  setting  the  proper 
sense  switch  to  the  "ON"  position.  Brieflv  the  fi\e  parts  are  concerned  with: 

(1)  printed  output  of  the  T(S),  TimeiV,  tnd  A(\)  arrays,  (:,  numerical  solution 
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Figure  3.  Flo**  Chart  for  Computation  of  Time  -  Temperature  Histories. 
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of  the  thermal  tissue  damage  integral  and  printed  output  of  results,  (3)  card 
punch  of  time -temperature  history  just  computed,  (4)  plot  of  time-temperature 
history  just  computed,  and  (5)  plot  of  time-temperature  history  for  surface 
temperatures  (x  *  o) .  Here  we  discuss  the  numerical  i  lution  of  the  thermal 
tissue  damage  equation.  The  others  are  discussed  under  output  in  the  Operating 
Instructions. 

Equations  3  and  4  yield  the  following  equation  for  tissue  damage 

rates  as  a  function  of  temperature  Tx, 

d n  ,  pe-AE/RTx  hq.  y 

at 

where  the  symbols  have  been  previously  defined.  The  values  ot  P,  Al. ,  and  R 
were  determined  as  follows  from  the  graph  in  Figure  4  (4); 

P  *  PI  •  2.1850  x  10-124 

and  AE/R  -  ER1  «  93,534,9  for  tissue  temperature,  Tx>  less  than  50°C,  and 

P  .  P2  »  1.823*  x  10*51 

and  Ab/R  -  ER2  *  39,109.8  for  tissue  temperature,  Tx>  equal  to  or  greater 
than  50*C . 

The  damage  rate  for  each  temperature  value  in  the  array  l(Ni  between 
TMINUP  and  TMINDN  is  computed  according  to  Eq .  9,  depending  on  whether  the  value 
of  the  temperature  is  less  than,  greater  than,  cr  equal  to  5o*C.  Values  of 
temperature  below  TMINUP  do  not  make  a  significant  contribution  to  the  total 
damage  Integra'.  The  heating  damage  integral  (for  values  of  t  <.  iture  between 

TMINUP  and  IPk),  HI,  and  the  cooling  damage  integral  .for  values  of  temperature 
between  TPk  and  TMINON),  Cl,  are  computed  according  to  the  trapetotdal  rule  for 
mtegrat  ion ; 
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Figure  4.  Daaaje  Rates  Derived  from  Radiative  Data 
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Eq.  10 


dt  dfl  .dO  .dO  -do  do 

r  3t,  *  23t2  *  23t3  .  *  23tE_i  *  ar. 


where  0  «  damage  integral 

dt  *  DT  «  time  interval  between  temperature  points 
dO/dtffi  *  damage  rate  at  given  temperature 
Following  the  computation  of  the  damage  integral  during  heating,  HI,  and  Cl, 
the  integral  during  cooling,  the  sum  HI  ♦Cl,  the  total  damage  is  stored  in 
FI.  Figure  S  is  a  flow  chart  of  the  integral  computation. 


OPERATING  INSTRUCTIONS 

General 

In  addition  to  a  deck  of  Hollerith  cards  containing  all  the  Fortran 
statements,  control  cards  are  necessary  for  the  operation  of  the  program.  How¬ 
ever,  since  the  number  and  format  of  the  control  cards  may  vary  somewhat  in  dif¬ 
ferent  computer  installations,  they  will  not  be  considered  here.  Details  on  the 
appropriate  control  cards  can  be  obtained  at  each  installation.  A  binary  deck 
containing  the  incremental  plotter  routine  to  operate  the  plotter  (7,8)  completes 
the  card  requirements. 

Various  modes  of  output  may  be  selected  by  means  of  sense  switches. 

The  operation  of  each  output  is  inhibited  by  placing  the  proper  sense  switch  in 
the  "ON"  position.  Figure  b  is  a  macro  flow  chart  of  the  logic  between  various 
sections . 

Input 

The  first  data  card  (Figure  7A)  contains  the  variable  L,  the  number  of 
graphs  to  be  plotted  or  the  number  of  t ime- temperature  histories  to  be  computed. 
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Figure  6 


Macro  Flow  Chart  of  Logic. 
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Figura  7.  Examples  of  Input  Data  Cards. 
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L  is  a  three-digit  integer  punched  in  columns  1-3  on  the  card  and  is  read  in  I 
format.  L  is  decremented  by  1  after  the  computation  of  each  history  and  tested 
for  zero.  When  L»0  the  program  is  terminated. 

The  second  data  card  (Figure  7ft)  contains  an  array  called  LILY  which 
has  12  four-digit  labels  used  to  label  the  ordinate  axis  of  the  graphed  outputs. 
The  first  value,  LILY(l),  is  punched  in  columns  1-4  on  the  card  in  the  form 
,35*  and  is  read  in  A  format  Successive  values  of  LILY  are  punched  in  successive 
columns  of  four,  thus  using  columns  1  through  48  for  the  entire  array. 

The  third  data  card  (Figure  7C)  contains  three  variables,  D,  B,  C,  used 
in  labeling  the  graphed  outputs.  They  are  the  date  expressed  in  month,  day,  and 
year  form,  thus,  00/00/00.  Each  value  is  three  digits  in  length  and  the  three 
variables  occupy  columns  1  through  9  with  the  following  form,  D  =  .00,  B=/00, 
and  C=/00,  The  values  are  read  in  A  format. 

The  fourth  data  card  (Figure  7D)  contains  the  values  PI,  P2,  ER1,  and 
ER2.  PI  and  P2  are  punched  in  columns  1-11  and  12-22  respectively  and  read  in 
E  format.  ER1  and  ER2  are  punched  in  columns  23-30  and  31-38  respectively  and 
are  read  in  F  format. 

The  remaining  data  cards  (Figure  7E)  contain  the  values  of  Q,  X,  DT, 

TAU,  TIMF(1),  Al,  A2,  ZER0TEMP,  end  PIESQRT.  All  the  values  are  read  in  F  format 
and  occupy  the  following  columns;  Q  in  columns  1-8,  X  in  columns  9-16,  DT  in 
columns  17-24t  TAU  in  columns  25-32,  TIME(l)  in  columns  33-40,  Al  in  columns 
41-49,  A2  in  columns  50-58,  ZEROTEMP  in  columns  59-66  and  PIESQRT  in  columns 
67-77.  Ail  values  cover  the  expected  range  of  the  values  with  a  sign  position. 

0u:put 

The  first  part  of  the  output,  the  printed  output  of  the  three  arrays 
T(N),  TIME(N),  and  A(N),  is  controlled  by  sense  switch  #1.  Figure  8  is  an 
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Figure  8.  Example  of  Array  Printout. 


example  of  the  array  printout.  The  values  of  Q,  Al,  A2,  and  TAU  are  printed 
at  the  top  of  the  page  for  identification  purposes,  followed  by  the  appro¬ 
priate  headings  for  each  cf  the  arrays.  The  arrays  are  each  split  in  half 
and  printed  in  six  columns  across  the  page.  It  will  be  noted  that  there  is  no 
variation  in  the  value  of  the  square  root  of  thermal  conductivity  during  the 
heating  phase  (0  <  TIME  <_  TAU),  while  the  square  root  of  thermal  conductivity 
during  the  cooling  phase  (TIME  >  TAU)  is  continually  decremented  by  0.0001  at 
each  time-temperature  point  as  mentioned  before. 

The  second  part  of  the  output,  the  printed  output  of  the  results  of 
the  numerical  solution  of  the  thermal  tissue  damage  integral,  is  controlled 
by  sense  switch  #2.  Figure  9  is  an  example  of  the  integration  printout.  The 
printout  consists  of  the  values  of  Q,  PI,  P2,  ER1,  ER2,  TPK,  TMINUP,  and 
TMINDN  for  checking  and  identification  purposes,  along  with  the  results  of 
the  integration  FI,  HI,  and  Cl.  It  is  seen  that  the  integration  printout 
directly  follows  the  array  printout.  Since  TPK  should  occur  at  TIME  =  TAU, 
TMINUP  should  be  the  value  of  temperature  nearest  to  44°C  during  the  heating 
phase,  and  TMINDN  should  be  the  last  value  in  the  T(x)  array,  these  values  can 
easily  be  checked  against  the  array  printout. 

The  third  part  of  the  output  punching  the  T(N)  array  on  cards  is  con¬ 
trolled  by  sense  switch  #3.  Figure  10  is  an  example  of  the  punched  data  cards. 
Each  array  is  preceded  by  a  card  containing  the  following  data:  PT  in  columns 
1-5,  NO  in  columns  6-8,  TMINUP  in  columns  9-15,  TMINDN  in  columns  16-22,  TPK 
in  columns  23-29,  and  Q  in  columns  30-35.  (Figure  10A)  The  array  from  TMINUP 
to  TMINDN  inclusive  is  punched,  eleven  values  of  temperature  per  card,  each 
value  having  a  maximum  of  six  digits  and  decimal  point,  three  places  to  the 
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*Jk 


U  1 

bdTMNur|M 

ra  am  iimiiii 


0 

0 


llllllll! 

‘immmmimQmmmiQimmmmmimiimmiim 

in 

IIIIIIIM 

•  inn 

II 1 1  It  1 1 

II  1  1  1  1  1  1  1  1  1  1  1  1 1 1 1 1 0 1 1 1 1 1 1 M 1 II 1  1  II  1  1  1  1  1  1  1  1  1  1  • 

1  1  1 

1 1 II II 1 1  1 

II  1  1  1  1 

1 1 II  tit  It  till!!! Ill  111  It  till  1  lllll  III llll  Ill  1 

in 

HIM!!!! 

;u  ii, 

iQllllll 

iiiQiinilt)  nniQiQniOOnnininn  i  innn 

in 

llllllll) 

nun 

iillilli 

iQlll  4  iQlllllllllll!  1  4  4  1  1  1  1  1  44  1  144444  4  III4I4 

4(4 

444441141 

4  4  4  4  4  4 

1 1 iQi 1 Qi 

IIIlQlltlllltt iQlllltlllllllllttlllllt  lllllll 

in 

nntitn 

lllll) 

llllllll 

IIIIIIIIMlIlQllllllllllltQllllllllllllllllll 

in 

iiiiiiin 

nun 

r  r 1 1 1 r  ii 

m  i  iQO" 1  >  in  i  iO>  >  >  i  >  >  >  >0* 1  n  i  ii  i  iii  i  it  1 1 1 1 1 

1  1  I 

1 1 1 1 1  •  1 1 1 

Mini 

IQllllll 

■iiQiiiiiiQiiiiiiQiiii'Qiiiiaiiiiiiiiiiiiiiii 

in 

•iiiiiiii 

•inn 

llllllll 

r  I  i  i  I  i 

Ilf  HttlltlfftttttllllltlMIIHf  till  tllttttlll 

i  <  •  j  •!  i  •  a  afr  r  i  *  i  »'•  »  i  •  i  o  n  a  '  ■  i  a  a  »•  •  f  ••  a*  j  i. 

in 

a  w  a 

iitimii 

i  jvwtktiMr 

mm 

4.  i  a  a  - 

*  w  \ 


0  0 
I  2  |  3 

nr  . 


0 

"  1 


0  0  0  0  0  1  0 

I  [2|3|4|5|6|7|el9|lO. 

i  n  nin  1 1 1 1 1 1  unn  n'n  n  1 1  i'i  1 1 1  m'n  n  1 1  i'm  m  1 1  •  i  •  ••  1 1 1 1 1  n  1 1 1 1 1  n  1 1 1 1 1  m  1 1 1 1 

iiiitiniOiiOiMiiiiiiiiMiiiiiiiimiii|inmiinmMiiiiigimiiiiiiiniii 

immiMiiimmMiiiMiIiiimtmignmtiOMiiiOMimmiiimntnOtDi 

nOO)nuiOniinOD>ii)iOr)iiiiOr>iiiiO>iO>>iD>>Ois>0<>>i>iOni>nO)i>>>iO>ni 

iOnnO»|MiniOHiiQiiiiiiiiiitiiiiiinininiiiiiiun!itninniOiinnnn 

mimmiimOiO'miiOmmiiimimhimmiimmiimihhiOmOmiiiiOiiiihOmim 

iiiOO|i"iO|,",iO-i*ii*0'OniOO,0|"OOn'iiiOiMiii|niiii|iiiiii|i  iOmiOiim 

,  ,4.1  nv  nnn  t  a  »»••»»  «  •  » a 4  4-  «  *  •  •  4  a i  v  »  «  «  a  m  a  t  a*  ii «  i  .  - 

t*  '  •  *  *  *  j  .  »  * 

B 


Figure  10.  Exaaple  of  Punched  Data  Cards. 
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right  of  the  decimal  point.  (Figure  10B) .  This  feature  was  originally  included 
in  the  program  to  provide  data  output  acceptable  for  use  in  another  program. 

The  fourth  part  of  the  output,  a  plot  of  the  entire  tissue  temperature 
array  T(N)  from  Tq  to  TMINDN  by  means  of  an  incremental  plotter  is  controlled 
by  sense  switch  #4.  Figure  11  is  an  example  of  a  tissue  temperature  plot.  The 
ordinate  of  the  plot  is  labeled  T(x)  and  has  a  range  from  25-75*C  with  the 
appropriate  labeling.  This  graph  procedure  is  used  to  plot  time-temperature 
histories  at  a  depth  x  below  the  surface  where  a  TPK  of  less  than  75*C  is 
des.rable.  Hence,  since  the  initial  surrounding  temperature  (TQ)  is  always 
32.5*C,  the  range  of  the  ordinate  is  sufficient  to  plot  tissue  temperatures  at 
depth.  The  length  of  the  abscissa  is  computed  for  each  plot  depending  on  the 
size  of  TIME(N).  If  TIME(N)  is  less  than  4.5  sec  the  length  of  the  abscissa 
is  5.0  sec  long;  if  TIME(N)  is  less  than  2.25  sec  the  length  of  the  abscissa  is 
2.5  sec  long.  If  TIME(N)  is  greater  than  4.5  sec,  the  length  of  the  abscissa 
is  equal  to  ten  times  the  following  truncated  value:  [ (TIME(N)/9. )♦! . ]  The 
abscissa  is  labeled  TIME (sec.)  and  covers  the  appropriate  range.  The  upper 
right  hand  comer  of  the  plot  contains  the  date  of  the  plot  and  the  values  of  Q, 
Al,  A2,  and  TAU  for  identification  purposes.  One  time-temperature  history  is 
plotted  per  each  graph. 

The  fifth  and  last  part  of  the  output,  a  plot  of  surface  time -temperature 
histories,  is  controlled  by  sense  switch  *5.  The  graph  is  exactly  the  same  as 
that  for  tissue  time-temperature  histories  except  that  the  range  of  the  ordi¬ 
nate  is  increased  to  25-90*C  to  handle  the  higher  TPK  of  the  surface  time- 

v 

temperature  histories. 

Since  there  are  five  sense  switches  each  with  "ON"  and  "OFF"  positions, 
there  are  thirty-two  different  modes  of  operation  of  the  program.  The  program 
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Figure  11.  Exaaple  of  Tissue  Teaperature  Plot. 


is  designed  to  operate  in  each  of  the  thirty-two  nodes;  however  only  a  few  of 
the  pssible  modes  are  used  in  actual  practice.  Some  of  these  are  briefly 
discussed  because  of  their  importance. 

1.  All  sense  switches  in  "OFF"  position  -  normal  mode  of  operation. 

The  three  arrays  T(N),  TIME(N),  and  A (N)  are  printed,  followed  by  a  printout 
of  the  results  of  the  tissue  damage  integral  as  shown  in  Figure  9.  The  T(N) 
array  is  punched  on  cards  according  to  the  format  shown  in  Figure  10  and  a 
graph  with  increased  ordinate  size  (2S*-90eC)  is  drawn  containing  both  surface 
temperatures  (the  plot  with  the  larger  TPK)  and  tissue  temperatures  (the  plot 
with  the  smaller  TPK)  as  shown  in  Figure  12.  bach  surface  time -temperature 
history  is  printed  out  as  shown  in  Figure  13  with  a  label  at  the  bottom  identify¬ 
ing  the  history  as  a  surface  temperature  history.  The  first  four  input  data 
cards  are  arranged  as  mentioned  before;  however  the  remaining  cards  containing 
values  of  Q,  X,  DT,  TAU,  TIME(l),  Al,  A2,  ZEROTEMP,  and  PIESQRT  are  arranged  in 
the  following  order.  Each  card  containing  values  of  Q,  X,  DT,  TAU,  TIME(l),  Al, 
A2,  ZEROTEMP  and  PIESQRT  for  a  time-temperature  history  at  some  depth  x  below 
the  surface  of  the  skin  is  immediately  followed  by  a  card  containing  exactly 

the  same  values  except  that  the  value  of  X  is  zero  (0.0).  Thus  the  value  of  L 
is  equal  to  the  number  of  such  paired  cards  or  the  number  of  graphs  drawn  but 
equal  to  one-half  the  number  of  time-temperature  histories  computed. 

2.  Sense  switch  »3  in  "ON"  position,  all  the  other  switches  in  "OFF" 
position  -  Operation  is  exactly  as  in  case  1  above  except  the  punching  of 
T(N)  array  is  inhibited. 

3.  Sense  switch  *5  in  "ON"  position,  all  the  other  switches  in  "OFF" 


position  -  Operation  is  exactly  as  in  Case  1  above  except  a  graph  with  normal 
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Figure  12.  Example  of  Surface  Temperature  and  Tissue  Temperature  Plot 
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Fifure  13.  Example  of  Surface  Time -Temperature  History  Printout. 


ordinate  size  (25*-7S*C)  is  drawn  containing  only  one  plot,  that  of  a  time- 
teaperature  history  at  depth  x  below  the  surface  of  the  skin.  The  value  of 
L  is  equal  to  maiber  of  data  cards  after  the  fourth  one  or  the  number  cf  graphs 
drawn  or  the  number  of  tiae- temperature  histories  computed. 

4.  Sense  switch  *4  in  "ON"  position,  all  other  switches  in  "OFF"  posi¬ 
tion.  Only  surface  tiae-teaperature  histories  are  computed  and  printed  out 

as  shown  in  Figure  13.  Integration  and  punching  of  surface  time-temperature 
histories  onto  data  cards  are  always  automatically  omitted  whenever  surface 
tiae-teaperature  histories  are  computed.  In  this  case  a  large  graph  is  drawn 
with  one  surface  tiae-teaperature  history  plotted  per  graph.  The  value  of  L 
is  determined  as  in  Case  *3. 

5.  All  sense  switches  except  *2  in  "ON"  position,  sense  switch  #2  in 
"OFF"  position.  The  only  operation  performed  is  the  evaluation  of  the  thermal 
tissue  daaage  integral,  the  results  printed  out  one  after  another  consecutively 
down  the  page. 

The  operational  analysis  of  the  other  nodes  of  operations  can  easily  be 
understood  by  reference  to  the  flow  chart  in  Figure  6.  It  should  be  noted  that 
the  computer  can  distinguish  between  surface  time-temperature  history  data 
(x-Q.O)  and  tine-temperature  history  at  depth  data  only  by  means  of  the  sequenc¬ 
ing  called  for  the  program.  Thus,  for  instance,  if  the  operator  ,oads  in 
surface  tiae-teaperature  history  data  and  places  sense  switch  »$  in  the  "ON"  posi¬ 
tion  the  computer  will  treat  the  data  as  tiae-teaperature  history  at  depth  data. 

It  is  the  responsibility  of  the  operator  to  make  the  input  data  consistent  with 
what  is  called  for  by  the  sense  switch  settings. 
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